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We investigate the high-pressure phase transitions in BiFe03 by single crystal and powder x-ray 
diffraction, as well as single crystal Raman spectroscopy. Six phase transitions are reported in the 
0-60 GPa range. At low pressures, up to 15 GPa, 4 transitions are evidenced at 4, 5, 7 and 11 GPa. 
In this range, the crystals display large unit cells and complex domain structures, which suggests 
a competition between complex tilt systems and possibly off-center cation displacements. The non 
polar Pnma phase remains stable over a large pressure range between 11 and 38 GPa, where the 
distortion (tilt angles) changes only little with pressure. The two high-pressure phase transitions at 
38 and 48 GPa are marked by the occurence of larger unit cells and an increase of the distorsion away 
from the cubic parent perovskite cell. We find no evidence for a cubic phase at high pressure, nor 
indications that the structure tends to become cubic. The previously reported insulator-to-metal 
transition at 50 GPa appears to be symmetry breaking. 
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I. INTRODUCTION 

Magnetoelectric multiferroics currently attract a great 
deal of interest, both because of intruiging coupling 
mechanisms between magnetism, ferroelectricity and fer- 
roelasticity, and because of the potential for new types of 
magnetoelectric devices [l]-[5j|. Among multiferroics, Bis- 
muth ferrite BiFeOs (BFO) is commonly considered as a 
model system [6|, and is perhaps the only material that 
is both magnetic and ferroelectric with a strong electric 
polarization at ambient conditions. Even though BFO 
attracts an important attention, the crystal structures of 
BFO as a function of temperature and pressure are still 
debated in the literature p]. 

Besides being a model multiferroic, BFO is also one of 
the very few perovskites that present both octahedra tilts 
@-[|[ and strong ferroelectric cation displacements, thus 
two structural instabilities, at room temperature. Ab- 
initio calculations suggest that in a number of perovskites 
both instabilities are present and compete fid Hi] . It is 
however important to remind that one instability usually 
prevails so that the vast majority of perovskites presents 
either tilts (SrTi03, LaA103 etc.) or ferroelectric cation 
displacements (BaTiOs, PbTi03 etc.) at a given tem- 
perature or pressure [8| , as schematized in Figure [U For 
some perovskites a crossover from one distortion to the 
other can occur as a function of an external parame- 
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ter, PbTiOs under high pressure is such an example fl2l |. 
BiFe03 is thus a rather exceptional case and turns out to 
be a model system for studying the competition between 
tilts and polar instabilities in perovskites. 

General rules for predicting phase transitions in per- 
ovskites are of longstanding interest. As a fruit of past 
studies it is now generally accepted that increasing tem- 
perature reduces both tilt instabilities (the tilt angle de- 
creases with increasing temperature) and polar instabili- 
ties (decrease of the polar cation displacement) (cl [l3l fl3 | . 
As a consequence, increasing temperature drives per- 
ovskites towards the parent cubic perovskite, although 
the ideal cubic structure might be in some cases beyond 
the decomposition or melting point. The effect of high- 
pressure on perovskites is more complex. Concerning oc- 
tahedra tilts, it was originally suggested by Samara et 
al. [15] that the phase transition temperatures T c of 
zone-boundary transitions in perovskites should always 
increase with pressure: dT c /dP > 0, i.e. the antiferro- 
distortive tilt angle increases with increasing pressure. 
However, ex per iments on LaA103 [16[ and later on other 
perovskites fl7H20| revealed that some perovskites de- 
crease their tilt angle and undergo phase transitions to 
higher-symmetry phases on increasing pressure. As a 
consequence, Samara's rule has been extended by An- 
gel et al. to a new general rule [l7[ which predicts the 
behaviour of octahedra tilts by taking into account the 
compressibilities of the different polyhedra, with so far 
no known exception. 

As to polar instabilities, the pioneering work by 
Samara et al. [15| describes that pressure reduces fer- 
roelectricity in ABO3 perovskites and even annihilates 




Figure 1: Schematic illustration of distortions and their com- 
bination in ABX3 perovskites. Each of the three inner circles 
represents one of the three main collective structural distor- 
tions in perovskites: (anti)polar cation displacements (blue), 
octahedra tilts (red) and first-order Jahn- Teller (green) dis- 
tortions, each illustrated by a schematic drawing and a mate- 
rial that adopts the distortion at ambient conditions. Materi- 
als outside the circles adopt the ideal cubic perovskite struc- 
ture exemplified by SrTiOs. BiFeOs is among the very few 
materials that adopts at 300 K both strong ferroelectric cation 
displacements and strong octahedra tilts. 



it for a critical pressure P c at which the crystal struc- 
ture becomes cubic. Early papers on BaTi03 [2l|, |22| , 
KNb0 3 0,0, PbTi0 3 [25|, |2H| amply confirmed this 
view and let to the widely accepted conclusion that po- 
lar perovskites adopt for P > P c a cubic Pm3m crystal 
structure. It was thus unexpected to find out in recent 
years that relaxor ferroelectrics do not stay/become cubic 
at high-pressure (e.g. Ref. [27-30]) and that ferroelectric 
PZT displays a wealth of phase transitions under high- 
pressure without becoming cubic either (e.g. Ref. |3l1 — 
l33j). However, the most striking observation was that the 
model ferroelectric and simple perovskite PbTi03 first 
reduces its ferroelectricity and becomes cubic but then, 
unlike commonly thought, becomes again non-cubic with 
tetragonal phases [HI, [3J, [35| through mechanisms in- 
volving oxygen octahedra tilting [l2[ and a theoretically 
predicted re-entrance of ferroelectricity (33436). 

The aim of our study is to go beyond the above dis- 
cussed perovskites with one instability only, by investi- 
gating the consequences of the competition of octahedra 
tilts and cation displacements on the crystal structure of 
BiFeOs, which is also magnetic. We focus on the effect 
of pressure, a parameter which allows driving to a very 
large extend structural instabilities, and thus their com- 
petition, thanks to otherwise unachievable large changes 
in cell volume and bond lengths. Despite literature re- 
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Figure 2: Schematic summary of the high-pressure investigs 
tions of BFO. 



ports on the effect of high-pressure on BFO, the actual 
symmetries in the low-pressure regime remain controver- 
sial and the very-high pressure regime has not been prop- 
erly investigated (see ref. |37l-[45| and following section). 
Here, we investigate BFO by using Raman spectroscopy 
and synchrotron X-ray diffraction on both powder and 
single crystals up to the very high-pressure of 55 GPa, 
a regime only rarely explored in functional oxides beside 
some notable exceptions fl2l l39l Efl |47| . 

We show that BFO undergoes a surprising richness of 
six phase transitions up to 55 GPa. Further to this, we re- 
port for BFO unusual large and puzzling unit cells under 
high pressure, which are unexpected for simple ABO3 
perovskites. The results are presented in light of the 
current literature on BFO and also in relation to other 
Bi-based perovskites. 



II. PRELIMINARY REMARKS 

The structure of BFO at high pressure has been investi- 
gated in the past in several studies, as schematized in fig- 
ure [2j Powder x-ray diffraction under high-pressure has 
been performed over the past few years by several groups 
with various outcomes. Gavriliuk et al. [H , EH and 
Zhu et al. [45| have reported no structural phase tran- 
sition in this pressure range. Haumont et al. [38] have 
identified two phase transitions at 3.5 and 10 GPa with 
a proposed phase sequence R3c — > C2/m — > Prima. 
Belik et al. |43[ have later confirmed the phase transi- 
tion at 4 GPa and identified an additional transition at 7 
GPa. They have identified the two intermediate phases 
as orthorhombic instead of the monoclinic C2/ra, with 
a transition sequence R3c — > Ortho I — > Ortho II — > 
Prima, and also pointed out reversibility issues. In a 
previous work, we have confirmed these transition pres- 
sures and emphasized the importance of good hydrostatic 
conditions j39[. Very recently, Kozlenko et al. [49] have 
reported a powder neutron diffraction study and identi- 
fied a high-pressure orthorhombic phase with space group 
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Exp. Sample 



PTM Beamline P range 



XRD 1 Single crystal Helium ID09A 0-12 GPa 

XRD 2 Single crystal Neon ID27 12-55 GPa 

XRD 3 Single crystal Helium ID27 30-55 GPa 

XRD 4 Powder Helium ID27 0-52 GPa 

Raman Single crystal alcohol 0-12 GPa 

Table I: List of the experiments reported in this paper with 
the pressure range investigated and the pressure transmitting 
media (PTM) used. 



Pbam between 3 and 8.6 GPa, distinct from the phases 
previously seen. In parallel, Raman spectroscopy has also 
been used to investigate phase transitions |37|, H] • Both 
studies revealed two phase transitions in the low-pressure 
range, but at somewhat different pressures, about 3 and 
9 GPa. Theoretical studies have also been devoted to 
the high-pressure behaviour of BFO (5(3453) and will be 
commented on in the discussion. 

Generally speaking, the determination of the struc- 
tures of BFO from single crystal diffraction in diamond- 
anvil cells (DAC) is complicated by (i) the limited access 
to reciprocal space, because of the geometry of the DACs, 
(ii) the complex domain structure that emerges as the 
crystal goes through the different phase transitions (iii) 
the very weak intensity of oxygen-related superstructure 
reflections (SSR) as compared to the heavier Bi and Fe. 
Powder diffraction on the other hand is not limited by 
the DAC geometry but overlapping bands, band broad- 
ening and a stronger background hamper the observation 
of weak reflections and small splittings. Both techniques 
can be used in a complementary way, although single 
crystal data give in general much better results, espe- 
cially at high pressure. In this paper, we report four 
different diffraction experiments carried out on powder 
and single crystal samples. The experiments are summa- 
rized in table U with the pressure ranges investigated. In 
the following, the diffraction results will be mostly taken 
from experiments 1 and 2 carried out with single crystals. 
Experiments 3 and 4 were carried out to check the repro- 
ducibility of the transitions; they will not be presented 
in detail. 



indicate that the crystals used are in a single domain 
state before the experiments. 

All experiments were performed in diamond-anvil cells 
(DAC). The diamonds have the Boehler-Almax design 
with cullets of 600 |Ltm (exp. 1 + Raman) or 250 ]im 
(exp. 2, 3 and 4). The pressure chamber was sealed by 
a stainless steel gasket. The pressure transmitting media 
(PTM) used are given in table HI 

X-ray diffraction experiments were performed on the 
ID27 and ID09A beamlines at the ESRF. At ID27, the 
beam was monochromatic with a wavelength of 0.3738 A 
selected by an iodine K-edge filter and focused to a beam 
size of about 3 |Ltm. At ID09A, the beam size is about 
20 [im and the wavelength (0.4144 A) was determined 
from the calibration using a standard silicon powder. The 
signal was collected in the rotating crystal geometry on 
a marCCD 345 (ID27) or mar 555 image plate (ID09A) 
with —30° < lj < 30° in 0.5 or 2° steps. A precise calibra- 
tion of the detector parameters was performed with a ref- 
erence silicon powder. The pressure was measured by the 
standard ruby fluorescence method [55|. The diffraction 
patterns from single crystal measurements were indexed 
with a home-made program based on the Fit2D software 
(56| . The refinement of the lattice constants from the 
peak positions was performed with the program UnitCell 
[57j . The CrysAlis software (Oxford Diffract ion- Agilent) 
was used for reconstruction of planes in reciprocal space. 
For the powder experiment, a gold powder was used as a 
pressure standard [58j . The pressure calculated from the 
gold standard was corrected for non-hydrostatic effects 
by following the procedure described by Singh [59] . 

The Raman spectra were recorded on a Jobin-Yvon 
Labram spectrometer with a low frequency cutoff at 100 
cm -1 . The exciting laser line was 633 nm. The laser 
power was kept at 10 mW on the DAC to avoid heating 
of the sample. The spectra are unpolarized. 



IV. RESULTS 



III. EXPERIMENTAL DETAILS 

The BFO powder was prepared by conventional solid- 
state reaction using high-purity (better than 99.9%) bis- 
muth oxide Bi2C>3 and iron oxide Fe203 as starting com- 
pounds. After mixing in stoichiometric proportions, pow- 
ders were calcined at T f = 820°C for 3 h. BFO single 
crystals were grown by the flux method. Details about 
the synthesis can be found in Ref. [54j. The single crys- 
tals used for the experiments were polished to a thickness 
of about 10 |Ltm with a lateral extension of 10 to 30 |Ltm. 
Both optical inspection under polarized light and XRD 



The analysis of the diffraction patterns reveals 6 phase 
transitions at 4, 5, 7, 11, 38 and 48 GPa with the 
phase sequence R3c — > Phase II — > Phase III — > 
Phase IV — > Prima — > Phase VI — > Phase VII. The 
R3c and Prima phases have been reported previously and 
are only confirmed by our experiment. In contrast, the 
structure of phases II, III, IV, VI and VII have not been 
determined. In the following, we divide the results into 
two distinct pressure ranges: the low pressure region from 
ambient pressure up to 15 GPa and the high pressure re- 
gion for the Prima phase and above. Finally we give 
the evolution of the lattice constants, volume and strains 
over the full pressure range. 



4 




5.5 GPa 

a = 5.484(1) A b= 16.674(1) A c = 7.737(1) A 




be 

































/ 3 



c = 2c c 



4.0 4.5 15 1 

2 theta 



17 18 100 200 300 400 500 

Raman shift (cm 1 ) 



Figure 3: (Left) Selected parts of the integrated diffraction 
patterns of the single crystal for the low-pressure experiment 
(A = 0.4144 A) with phases R3c (green), II (blue), III (or- 
ange), IV (black) and Pnma (red). (Right) Raman spectra 
up to 12.1 GPa. 



Figure 4: Maps of (hkO), (2kl) and (hOl) reciprocal planes 
obtained at 5.5 GPa. Due to the orientation of the sample 
in the DAC, we have almost no access to the (Okl) plane 
and the sampling for (2kl) is less precise. The relations of 
the orthorhombic unit cell axes to the simple cubic axes are 
recalled. 



A. Low pressure region: 0—15 GPa 

An overview of the experimental evidence for the 4 
phase transitions in the low pressure range is given in 
figure [3] where we present the integrated single crystal 
diffraction pattern (left) and the Raman spectra (right). 

The transition from the low-pressure R3c phase to 
phase II is marked by the emergence of new superstruc- 
ture reflections (blue in figure [3]). In addition, splittings 
of diffraction peaks show the emergence of a complex 
domain structure. A simple observation under a micro- 
scope also shows that the single domain has split into a 
multitude of very small domains in phase II. 

At 5 GPa, we observe the disappearance of many 
diffraction peaks, revealing a simplification of the domain 
structure, and a marked change of the Raman signature, 
marking the transition from phase II to phase III (or- 
ange in figure [3]). The transition from phase III to phase 
IV (black in figure [3j) at 7 GPa is not characterized by 
pronounced change in the diffraction pattern, in terms 
of emergence or disappearance of diffraction peaks. A 
close inspection of the Raman spectra does reveal small 
changes at 7 GPa, notably at low wavenumbers with the 
emergence of a mode at 100 cm -1 . The most clear sig- 
nature of this transition are sudden shifts of diffraction 
peaks at low 26 toward lower angles, which can only be 
associated with sudden changes in the lattice constants. 
This change coincides well with the transition seen at the 
same pressure in the powder diffraction experiment by 
Belik et al. [43|. At 11 GPa, we observe the transition 
to the Pnma phase by changes in both the diffraction 
patterns and the Raman signature (red in figure [3]) . 



The diffraction patterns for all three phases can be in- 
dexed with an orthorhombic unit cell, taking into account 
the complex domain structure. This cell has lattice con- 
tants a Q = 5.484(1) A, b Q = 16.674(1) A, c Q = 7.737(1) A 
at 5.5 GPa and a comparatively large volume (Z = 12). 
For this cell, pseudo-cubic lattice constants can be cal- 
culated as a pc = a Q /\/2, b pc = b Q /(3\/2) and c pc = c Q /2. 
In figure |4j we reproduce maps of reciprocal space at 
5.5 GPa and give the relations between the lattice vec- 
tors of the orthorhombic cell and those of the parent cubic 
perovskite. 

Metrically, this cell is close to the orthorhombic cell 
proposed by Belik et al. [43], except that we find a c 
axis doubled in all phases, as demonstrated at 5.5 GPa 
by the presence of the (2fcl) reflections in figure HI The 
pattern can also be indexed with the C2/m phase pro- 
posed by Haumont et al. [38], but the diffraction does 
not give evidence for a monoclinic distortion. On the 
other hand, the orthorhombic cell identified by Kozlenko 
et al. in their neutron diffraction study UK , with lattice 
constants V2a pc x 2y/2a pc x 2a pc , cannot explain the ob- 
served patterns. 

The inspection of the reflections in figure [4] shows the 
typical conditions for / centering. In addition, the inspec- 
tion of the (hOl) plane shows that only reflections with 
h,l = 2n are present. This leaves the possibility for two 
extinction symbols / — (ac)— and / — c6, which we can- 
not distinguish here because of an insufficient coverage 
of the hkO plane due to the DAC geometry. This leads 
to 6 different possibilities among 4 space groups: Ima2 
(46), I2cm (46), I2cb (45), Imam (74), Imcm (74), and 
Imcb (72). 
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Within this choice of space group assignments, we can 
tentatively infer information from the evolution of the 
domain structure at the transition II — >III. For a ferroic 
transition between two orthorhombic space groups, the 
formation of twins can only occur for a transition from 
point group mmm for the parent phase to mm2 or 222 
|60| . The simplification of the domain structure therefore 
suggests that phase II has a space group Jma2, 12cm or 
I2cb while phase III has Imam, I mem or Imcb. For the 
transition III — )TV, on the other hand, the diffraction 
patterns does not show any sign of a symmetry change, 
nor any change in the domain structure. This leaves the 
possibilities of a transition between any of the three cen- 
trosymmetric space groups mentioned above, which in- 
cludes isostructural transitions. 



B. High-pressure region: 11—60 GPa 

In the high-pressure region, we report two phase transi- 
tions at 38 and 48 GPa. Each phase will be separately de- 
scribed. Finally, we discuss the issue of non-hydrostatic 
stress conditions. 



1. Phase V: 11-38 GPa (Pnma) 

Between 11 and 38 GPa, we find the now well es- 
tablished GdFe03-type Pnma phase. A representative 
diffraction pattern is shown in figure [5j In this phase, 
the crystal shows peak splitting associated with a four- 
variant domain structure, as is commonly observed for 
Pnma crystals. As the pressure increases however, the 
intensity ratio changes and reflects the formation of a 
quasi single domain state (figure [5] (b) and (c)). 
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Figure 5: (a) Diffraction pattern of a BiFeOs single crystal in 
the Pnma phase integrated over the full angular range. The 
peaks show a splitting associated to the domain structure, 
(b) Evolution of the (222) c reflection under pressure and (c) 
selected corresponding 20 scans. 



2. Phase VP. 38-48 GPa 

At 38 GPa, additional sets of superstructure appear as 
can be seen in figure [6] on the reciprocal space map of the 
(hkO) plane. Metrically the new cell can be derived from 
the Pnma cell by a doubling of the a-axis. At 43 GPa, 
we find a Q = 10.546(4) A, b Q = 7.393(3) A and c Q = 
5.113(16) A. For this phase, we have Z = 8 and we cal- 
culate pseudo-cubic lattice constants as a pc = a Q /(2\/2), 
frpc = b Q /2 and c pc = c Q /y/2. Note that diffraction spots 
of the (100) reflexions are clearly visible, but (200) are 
not seen, which can only be explained by their weak in- 
tensity. Also, in this pattern, a small portion of a second 
crystal is visible, but is disoriented with respect to the 
main crystal, which may indicate that it comes from a 
piece of the crystal broken rather than a phase transition 
induced twin. The intensity of this additional peaks is 
however much weaker than the Bragg peaks of the main 
crystal and makes any confusion unlikely. 

The analysis of the reflection conditions leads to two 
possible extinction symbols (Pb and Pn ). These 



two possibilities differ by a condition 0kl:k and 0kl:k + I 
respectively. Unfortunately, our collected single crystal 
pattern does not allow to distinguish the two possibili- 
ties because of the restricted coverage of the Okl plane 
due to the experimental geometry and the orientation of 
the crystal in the cell. We found that the weak (100) 
reflections, although unambiguous in the single crystal 
pattern, cannot be seen in the powder diffraction pat- 
tern. The ratio of the intensities of these weak super- 
structure to the main Bragg peaks is about 1/10000. 
Powder diffraction alone would fail to identify the cor- 
rect unit cell and did not allow to discriminate between 

the two possibilities Pb and Pn . We are left 

with 5 possible space groups: the 3 polar groups Pb2\m 
(26), Pbm2 (28) and Pnm2\ (31) and 2 non polar groups 
Pbmm (51) and Pnmm (59). Based on the generally ac- 
cepted fact that ferroelectricity vanishes under pressure 
and that there is no known ferroelectric perovskite be- 
yond 15 GPa, we consider the last two possibilities as 
the most probable. In addition, we note that the Pnmm 
opens the possibility for a group-subgroup relation with 
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Figure 6: Reciprocal space maps of the (hkO) plane at 30 GPa 
(Phase Prima), 43 GPa (Phase Pnmm) and 53 GPa (Phase 
Cmcm), taken from experiment 2 (Neon as PTM). The im- 
ages have been rescaled for convenience so that the overall 
volume reduction is not apparent. The lattice constants and 
reciprocal lattice vectors for the different cells are recalled, 
including the monoclinic axes for the very last cell. 



the Prima phase. 



3. Phase VII: P > 48 GPa 

The transition from phase VI to phase VII is evidenced 
by a clear change in the diffraction pattern. A recipro- 
cal space map is presented in figure [6] This pattern can 
be indexed in a C-centered orthorhombic cell with the 
same axis directions as in phase VI. This new orthorhom- 
bic cell is derived from the Pnmm cell by a doubling of 
the b axis leading to lattice constants a = 10.417(4) A, 
b = 14.516(6) A and c = 4.961(15) A at 53 GPa, with 
16 formula unit per unit cell. With this indexing, the 
analysis of the reflection conditions leads to the extinc- 
tion symbol C—c— which allows for three different space 
groups: Cmc2i (36), C2cm (40) and Cmcm (63). The 
first two space groups are polar. With the same argu- 
ment as before, we regard Cmcm as the most probable 
space group for this phase. 

Alternatively, the pattern can be indexed in a mono- 
clinic cell. This cell has a volume half of the orthorhom- 



bic cell (Z = 8) and reciprocal lattice vectors related to 
the orthorhombic by a^ = + b£>, = a£> - b^, 
= c o> as depicted in figure [6] The orthorhombic and 
monoclinic cells are equivalent in the limit a^ = c* M . 
The analysis of the extinctions leave the possibilities for 
5 monoclinic space groups for which no reflection con- 
dition applies: P2 (3), P2 X (4), Pm (6), P2/m (10), 
P2\jm (11). Discarding as before the polar space groups, 
the two space groups P2/m and P2\/m appear as more 
likely. 



4- Importance of non-hydrostatic stress 

As long as the pressure-transmitting medium (PTM) 
remains liquid (P < 12 GPa for He, 4 GPa for Ne, 
10 GPa for the classical 4:1 methanol-ethanol mixture), 
the stress field can be safely regarded as hydrostatic. 
However, above this limit, a so-called deviatoric compo- 
nent may add to the hydrostatic stress field. The devia- 
toric stress has a rotational symmetry around the loading 
axis of the cell, with a stronger compressive stress along 
the axis, as checked by Zhao et al. [61] in an experiment 
with silicon oil as PTM. Non hydrostatic stress should 
then be regarded as unavoidable at high enough pressure, 
even though they are very often (legitimately) neglected 
when noble gases, particularly helium, are used as PTM. 
It is nonetheless of concern, for a meaningful comparison 
between different experimental works and with ab-initio 
calculations, to consider how such non- hydrostatic stress 
affects the outcome of an experiment. This can be done 
for example by comparing experiments performed in dif- 
ferent PTM or, for single crystal investigations, by com- 
paring experiments with crystals of different orientations. 

In a previous diffraction work on BiFeOs single crys- 
tals in the 0-10 GPa range, we have shown that strongly 
non-hydrostatic stress have a remarkable effect on the 
phase transitions, by stabilizing a structure that is not 
observed in hydrostatic conditions [39|. Here, we have 
shown an evolution of the domain structure in the Prima 
phase (figure [5]), which can be explained by the onset 
of non-hydrostatic stress favoring one particular domain 
orientation. However, we also make the following obser- 
vations. First, the same phases were observed in two ex- 
periments performed with different pressure-transmitting 
media (helium vs. neon). Moreoever, the powder diffrac- 
tion experiment shows the same transitions, although 
with sligthly different pressures. Therefore, we believe 
that non-hydrostaticity of the stress field plays a minor 
role and is only likely to shift transition pressure and af- 
fect slightly the evolution of the lattice constants, but 
without affecting the structure of the phases. 

In addition, complications may arise from the ferroe- 
lastic character of BFO. As the crystal undergoes dif- 
ferent phase transitions, ferroelastic domain formation 
and reorientation may act as stress relief and modify the 
stress state of the sample. This can make hydrostatic 
conditions difficult to be achieved, even at low pressures. 
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Figure 7: Evolutions of the pseudo-cubic lattice constants 
(top), volume (middle) and total strain etot (bottom) in the 
full pressure range investigated. Data are from the single 
crystal experiments 1 and 2. In the volume plot, the data 
from the powder diffraction experiment have been added for 
comparison, and the volume predicted by the equation of state 
has been extrapolated over the full pressure range as a guide 
to the eye. 



Ferroelasticity in BFO is unfortunately not extensively 
studied in literature, but deserves further attention in 
the light of the above comments. 



C. Volume, lattice constants and strain 

In figure we present the evolution of the lattice con- 
tants, volume and strain over the full pressure range. 
The strain is calculated as follows. As a reference state 
for all the phases, we choose a cubic cell with lattice 



constant ao = (V/Z) 1 / 3 . The general expressions for 
the strain components can then be found in (r32| . 
In the orthorhombic system, the expressions reduce to 



1, e 2 2 = b pc /a - 1, e 33 



/ao - 1, 



^pc? ^pc and c pc are the pseudo-cubic lattice con- 



en — ttpc/^0 

where a pc , b v 

stants. We then calculate the "total" strain using Aizu's 
expression e to t = \l^2 e %- This choice is somewhat arbi- 



trary, but enables us to quantify and compare the strain 
in the different phases. 

In the evolution of the lattice constants, the most 
remarkable feature is the transition III — >YV that is 
marked by a "switching" of the a pc and c pc parame- 
ters for a volume (and a strain) that shows no significant 
jump, as already reported by Belik et al in their powder 
diffraction study [43|. We have no precise decription of 
the microscopic mechanism underlying this phase tran- 
sition, but we note that a behavior of lattice constants 
similar to the present case was observed in BiMnOs [63[ 
under pressure, where it was associated to a change in 
the magnetic order. 

The evolution of the strain reveals the different distor- 
tion of the successive phases. With increasing pressure, 
the structure tends first toward a metrically cubic struc- 
ture, whereby the strain approaches zero. This evolu- 
tion is in good agreement with the literature [38]. The 
strain does not change appreciably through the transi- 
tions II — ?TII — )TV. Then, it shows a jump at the IV — > 
Pnma, which appears more distorted than the original 
R3c structure, and remains constant within experimental 
uncertainty through the Pnma stability range. At higher 
pressures, above 40 GPa, the tendency of an increase of 
the distortion away from a cubic metricity is very clear. 

In the volume plot, we have also plotted the values de- 
termined from the powder XRD experiment. The volume 
measured from the powder appears to deviate slightly 
from the values measured on the single crystal, with 
might be due to the different pressure-transmitting me- 
dia used (helium vs. neon). A fit of the pseudo-cubic 
volume in the stability region of the Pnma phase with a 
2 nd -order Birch-Murnaghan equation leads to Vi2GPa = 
56.43(8) A 3 and K 12GPai = 218(5) GPa for K[ 2GPa being 
fixed to 4. This value compares reasonably well with the 
previously published values of 238 GPa [38] , 

Last, we discuss the order of the transitions. From 
the discontinuities in the volume and lattice constants, 
it is clear that the transitions R3c — >II, III— )>IV and 
IV— » Pnma at 4, 7 and 11 GPa, respectively, are of 
first order. Transitions II— >-III and Pnma — >• Pnmm 
at 5 and 38 GPa, on the other hand, present possible 
group-subgroup relationships and show no appreciable 
volume discontinuities, rather suggesting second-order 
transitions. This also seems to be the case for the last 
transition Pnmm — >• Cmcm transition at 48 GPa, al- 
though it is difficult to be conclusive, due to the un- 
certainty on the volume measurements and the limited 
number of points close to the transition. 
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V. DISCUSSION 
A. Multiple phases of BFO 

Perovskites are known for their richness of structural 
and physical phase transitions with different possible in- 
stabilities which can coexist and compete. Both, a funda- 
mental and technological interest derive from this charac- 
teristic. One consequence of the various instabilities is an 
often complex phase diagram of perovskites, for instance 
in the P-T or x—T space which are the most explored dia- 
grams. Besides perovskite-type manganites, ferroelectric 
perovskites are among the most investigated materials 
with the representative examples PbZri-^Ti^Os (PZT) 
for the x-T space (gi^ and BaTi0 3 (BTO) or KNb0 3 
(KNO) for the P-T space @,[67j]. Both examples display 
complex phase diagrams. Nevertheless, it is important 
to realize that when only one parameter (P, T or x) is 
changed no more than 3 phase transitions are observed, 
e.g. the rhombohedral-orthorhombic-tetragonal-to-cubic 
phase sequence in BTO or KNO under pressure (24l |67|. 
Most other materials display fewer phase transitions un- 
der pressure: e.g. CaTiOs (no transition [47]), SrTiOs 
(one transition |4fj), or PbTi03 (two transitions [12]). 
In the light of this, the here observed 6 phase transitions 
of BiFeOs are remarkable. It is plausible to relate this 
observation to the fact that BFO presents at ambient 
conditions the rare coexistence of octahedra tilts, cation 
displacements and magnetic order. 

A recent study by Dieguez et al. (52| reinforces this in- 
tuitively expected intrinsic richness of structural phases 
in BFO by performing a systematic search for poten- 
tially stable phases by using first-principles methods. 
As a major outcome of their work, it has been demon- 
strated that BFO can present an unusual large num- 
ber of (meta) stable structures of which the balance or 
stability is determined not only by the traditional soft 
modes but also by secondary modes. Our observation 
of a pressure-induced rich phase sequence corroborates 
experimentally the predicted structural richness of BFO. 
However, our experimental results go beyond the predic- 
tions in the sense that we report new phases with unusual 
large cell dimensions, illustrating an even more complex 
energ y la ndscape than suggested by Dieguez's calcula- 
tions [52[ which were restricted for computational reasons 
to smaller cells. The physical origin of such large unit 
cells is not yet understood, although its systematic occur- 
rence in the low-pressure regime of Bi-based perovskites 
such as BiMnOs, BiScOs and BFO has been hypothet- 
ically related to the presence of the so-called lone-pair 
in bismuth (see discussion in Ref. (HI). We note how- 
ever, that we have shown here that BFO presents such 
large cells also at very high pressure, in sharp contrast 
to BiMnOs where the unit cells at very high-pressure 
are small [68j |. This suggests that the mechanism for 
very large unit cells has to be understood beyond the 
simple presence of Bi and thus deserve further atten- 
tion, namely through adapted first-principle calculations. 



NaNb03 [69] and AgNb03 (z3| show similarly large unit 
cells that have been interpreted in terms of complex tilt 
systems. Here, we understand by "complex tilt systems" 
collective tilts of the oxygen octahedra, whose description 
requires more than a doubling of the lattice constants of 
the cubic primitive perovskite cell and as such go beyond 
the now classical approach by Glazer [7, 71]. These com- 
pounds also show a coexistence between octahedra tilts 
and (anti) ferroelectric cation ordering and might provide 
some understanding. 



B. Stability and compressibility of the Prima phase 

A particular phase within this complex structural land- 
scape is the orthorhombic Prima phase. Generally speak- 
ing, Prima is one of the very common structures among 
perovskites and is known to withstand very high pres- 
sures. In particular, the rare-earth orthoferrites REFeOs 
maintain their Prima structure at high-pressure, even 
through their large volume collapse associated with a 
spin transition around 50 GPa [72|. In BFO, it presents 
a remarkable stability from 11 GPa up to 38 GPa. This 
constrasts with the effect of high temperature where it 
persists only in a narrow temperature range (73| . The 
stability is in agreement with calculations [52[ that have 
predicted that the Prima phase constitutes BFO's most 
stable phase beside the room temperature R3c phase. 

In the Prima phase, the evolution of the tilt angles 
under pressure is of particular interest due to the impor- 
tance of these angles for the electronic band structure. 
The two tilt angles Ti-Ol-Ti (a ± ) and Ti-02-Ti (a 2 ) that 
characterize the Prima structure [8| are best determined 
from the atomic positions obtained by Rietveld refine- 
ments, but the quality of our powder data does not allow 
to determine atomic position with satisfying uncertain- 
ties. Tilt angles can also be calculated from the lattice 
constants under the assumption of undistorted polyhe- 
dra [8|, but these formulas are known to give only rough 
estimations of the tilt angles [47|. From these formula 
and our measured lattice constants, we can estimate that 
the octahedra tilt angles decrease in the Prima structure 
by less than 2° between 10 to 38 GPa, with only little 
changes in the spontaneous deformation (figure [7j) 

Following Angel et al. compressibility rules [17( the tilt 
angle is predicted to decrease with increasing pressure 
for most A 3+ B 3 ^~Os perovskites, while it increases for 
A 2+ £> 4+ 03 perovskites. Our experimental observation 
appears in qualitative agreement with this rule, in spite 
of the presence of the lone pair of bismuth. Such small 
changes have been reported [19| to occur in orthoferrites 
in the presence of large and heavy A cations, as is the case 
here for BFO. In future works, it would be interesting to 
investigate in more detail rare-earth orthoferrites in a 
similar manner to obtain a better understanding of the 
particular role of bismuth on the compressibility of the 
perovskite. 
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Figure 8: Revised schematic pressure-temperature phase di- 
agram of BFO. The symbols are experimental points taken 
from this work (open symbols) and Refs. 0, [?3| (full sym- 
bols). The transition temperature to a cubic phase at ambi- 
ent pressure lies beyond the decomposition temperature and 
is a theoretical prediction [TtJ . 



C. Implications for the insulator-to-metal 
transition 

Another aspect of the very high-pressure regime is the 
reported occurrence of an insulator-to-metal (MI) phase 
transition accompanied by the disappearance of magnetic 
order around 50 GPa j40M42Ll48| |. There is an ongoing dis- 
cussion in the literature on the origin of this MI transition 
in terms of a high-spin to low-spin crossover, a Mott-type 
transition and related changes in crystal symmetry or not 
d, [H, l53l. I71 l75j. Our structural study does naturally 
not allow discussing in detail the proposed models of this 
physical phase transitions. However, our study provides 
evidence that the MI transition is symmetry breaking, 
and it rules out a transition to the earlier considered 
cubic [Hi or rhombohedral [H, [53[ phases. Rare-earth 
orthoferrites REFeOs present similar transition around 
50 GPa which has been assigned to a low-to-high spin 
phase transition of the Fe 3+ cations accompanied by a 
steep anomaly in the cell volume, particularly marked in 
LuFeC>3 f72ll76l|. Such a steep anomaly, characteristic for 
low-to-high spin phase transition is not observed in our 
measurements, thus rather questioning such a scenario. 



VI. CONCLUDING REMARKS 

In this study we have presented a structural investiga- 
tion of the model multiferroic perovskite BiFeOa under 
high pressure up to 60 GPa. The complementary use of 
synchrotron X-ray diffraction and Raman scattering re- 
veals a remarkable richness of 6 phase transitions in the 
investigated pressure range. The occurrence of these six 



pressure-induced phase transitions together with the four 
transitions as a function of temperature reported in the 
literature d, [HH, [ll, |7§| suggests a very complex phase 
P-T phase diagram that is characterized by the compe- 
tition of the BFO-characteristic instabilities: octahedra 
tilts, cation displacements, magnetic order and MI phase 
transitions. 

In a first attempt, Catalan and Scott [6] have sug- 
gested in 2009 a preliminary and very useful phase di- 
agram. Since then, our understanding of the different 
phase transitions has increased and we propose in Fig. 
[8] an updated schematic phase diagram, which is based 
on both literature results and the outcome of our present 
pressure work. In the following, we will discuss the dif- 
ferent features of this phase diagram with a particular 
emphasis on the essential points of our present work. 

(i) In agreement with the calculations by Dieguez et 
al. [52], the phase diagram is characterized by two 
particularly stable phases. First, the rhombohedral 
R3c phase (orange in Fig. [8]) which is at ambient 
pressure stable over 1000 K. Despite this temper- 
ature stability, the R3c structure is rapidly desta- 
bilized under hydrostatic pressure and disappears 
at a modest pressure of 4 GPa. Second, the or- 
thorhombic Pnma structure (blue), which is stable 
at 300 K between 11 and 39 GPa. We have shown 
that the Pnma phase is characterized by only small 
changes in the spontaneous deformation (and octa- 
hedra tilts), so that the compression mainly acts 
through bond shortening. 

(ii) Between the stable R3c and Pnma phases, we ob- 
serve under high-pressure three different phases of 
orthorhombic symmetry, thus a rather complex en- 
ergy landscape within only a few GPa. All three 
intermediate phases are characterized by unusual 
large unit cells, implying complex tilt systems and 
possibly off-center cation displacements triggered 
by the lone pair of bismuth. In the phase dia- 
gram we consider this intermediate region (striped 
in Fig. [Hj) as only one single region that is charac- 
terized by a multitude of competing phases, which 
might well reveal even more richness when varying 
the parameter temperature. The suggestion of a 
distinct antipolar phase by Kozlenko et al. [49[ 
supports this view. We also remind that new 
phases can easily be induced by non-hydrostatic 
conditions in this region [39| . The phase diagram is 
consistent with the idea of a complex deformation- 
temperature phase diagram as can be also explored 
by temperature-dependent investigation of strained 
thin films (although the biaxial strain in thin films 
can of course not be directly compared to hydro- 
static conditions). The very recent reports [80-83] 
of a structural phase transition close to room tem- 
perature in highly compressively strained BFO thin 
films adds further support to this. 
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(iii) The very high-pressure regime is characterised by 
two further phase transitions at 38 and 48 GPa, 
which are again characterized by large unit cells. 
Interestingly, we observe in this regime an increase 
of the distortion away from the cubic parent per- 
ovskite cell, which contrasts the earlier sketched (6| 
high-pressure tendency towards a cubic structure. 
On the other hand, this is consistent with ab-initio 
calculations by Gonzalez- Vazquez and Ihiguez ac- 
cording to which the simple cubic phase does not 
become the ground state even at high pressure [53] . 
At high temperature, the most recent studies [73] 
show that the cubic phase can only be expected at 
temperatures higher than the decomposition point 
and is therefore regarded as unstable. 

(iv) The insulator-to-metal phase transition in BFO is 
symmetry breaking. However, contrary to the ear- 
lier phase diagram @ the pressure-induced param- 
agnetic metallic phase is not cubic, but rather of 
orthorhombic symmetry. 



The above discussion and the phase diagram in Fig. [8] 
is based only on temperature measurement at 1 bar and 
pressure measurements at 300 K. As a consequence, the 
unexplored large intermediate P-T region remains hy- 
pothetic and is likely to give rise to new and unexpected 
features. We strongly encourage the exploration of the 
strain-temperature phase diagram, be it by hydrostatic 
pressure or strain in thin films, and particularly into the 
mid-pressure region with competing phases. 
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